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Abstract — A multi-access multiple-input multiple-output 
(MIMO) relay differential chaos shift keying cooperative 
diversity (DCSK-CD) system is proposed in this paper as a 
comprehensive cooperation scheme, in which the relay and 
destination both employ multiple antennas to strengthen the 
robustness against signal fading in a wireless network. It is 
shown that, with spatial diversity gains, the bit error rate (BER) 
performance of the proposed system is remarkably better than 
the conventional DCSK non-cooperation (DCSK-NC) and DCSK 
cooperative communication (DCSK-CC) systems. Moreover, 
the exact BER and close-form expressions of the proposed 
system are derived over Nakagami fading channels through 
the moment generating function (MGF), which is shown to be 
highly consistent with the simulation results. Meanwhile, this 
paper illustrates a trade-off between the performance and the 
complexity, and provides a threshold for the number of relay 
antennas keeping the user consumed energy constant. Due to 
the above-mentioned advantages, the proposed system stands 
out as a good candidate or alternative for energy-constrained 
wireless communications based on chaotic modulation, especially 
for low-power and low-cost wireless personal area networks 
(WPANs). 

Index Terms — multiple-input multiple-output (MIMO), relay, 
bit error rate (BER), differential chaos shift keying cooperative 
diversity (DCSK-CD), Nakagami fading channel. 



I. Introduction 

Inspired by the wideband nature, chaotic signals have shown 
a great promising for spread-spectrum communication systems 
0J-0. Among the digital communication schemes proposed 
so far, the chaos shift keying (CSK) and differential CSK 
(DCSK) have been widely studied and analyzed in recent years 
EJ— El attracted by their excellent capacity against multipath 
fading and time delay. Between these two schemes, DCSK is a 
non-coherent method which only requires frame or symbol rate 
sampling instead of channel estimation J9J. Furthermore, the 
hardware advantage of the DCSK scheme does not demand 
extra spreading and dispreading circuits in the mono-user 
case 1101 . As a result, this method is more practical at the 
expense of some performance degrading in comparison with 
its coherent counterpart. 
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In a wireless network, signal fading aroused from multipath 
propagation is one important factor that influences the system 
performance. Many methods have been proposed to overcome 
this weakness. On the one hand, some chaotic cooperative 
communication systems have been suggested to deal with 
multipath fading, e.g., the DCSK cooperative communication 
(DCSK-CC) systems ATI . Ifl2l . which can provide a consider- 
able gain over the conventional DCSK systems. On the other 
hand, multi-antenna techniques are particularly attractive as 
they can be easily combined with other forms of diversity. 
The use of multiple antennas at all terminals or only the 
receiver ends of a wireless link, which are called multiple- 
input multiple-output (MIMO) or single-input multiple-output 
(SIMO) techniques, demonstrates a significant improvement 
in link reliability and spectral efficiency iTOl . lfT4l . Motivated 
by these observations of advantages, several MIMO lfl5l . Ifl6l 
and SIMO systems 11171 have been carefully studied aiming to 
improve the performance of the DCSK communications over 
fading channels. As a potential application, applying MIMO 
technique to relay networks has been considered in recent 
years lfl8l . |[T9l , referred to as cooperative diversity. Moreover, 
some theoretical research on MIMO relay channels has been 
carried out, such as their capacity analysis lEUl . 11211 . 

Due to the advantages of low cost and low complexity, 
some ultra-wideband (UWB) systems based on DCSK or 
FM-DCSK have been considered for wireless personal area 
networks (WPANs) and wireless local networks (WLANs) 
G2-||24). In addition, a SIMO FM-DCSK UWB system has 
been proposed (25] to further improve the system performance 
in such environments. 

However, little research has been conducted focusing on 
the performance improvement of the multi-access DCSK sys- 
tems. To address this problem, the MIMO relay technique is 
introduced to the multi-access DCSK system in the present 
paper, which will be called an MIMO relay DCSK-cooperation 
diversity (DCSK-CD) system. Although channels are assumed 
to be Rayleigh fading in most studies for DCSK systems, it 
was pointed out that path fading statistics are more suitably 
to be described by a Nakagami distribution l26l . In fact, 
Nakagami channels have already been used for research on 
cooperation systems 11271 . ll28l . Motivated by this observation, 
the presently proposed system will be investigated over fading 
channels subjected to Nakagami distribution. In the proposed 
system, the relay and destination can employ multiple antennas 
to construct two SIMO and one MIMO systems of DCSK, 
on the source-relay link, source-destination link and relay- 
destination link, respectively. The new cooperation scheme can 
significantly enhance the capacity against signal fading in a 
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wireless multipath fading channel with moderate complexity. 
The performance of the proposed system is much better 
than the conventional DCSK-NC and DCSK-CC ones fl2) 
according to our theoretical analyses and simulated results, to 
be reported below, where more gain can be achieved as the 
number of antennas used increases. Moreover, a threshold for 
the number of relay's antennas will be derived provided that 
each user is subjected to an energy constraint. Furthermore, 
a close-form BER expression over Nakagami fading channels 
will be derived by using moment generating function (MGF), 
which agrees well with the numerical simulations. 

The remainder of this paper is organized as follows. In 
Section HJ the principles of DCSK modulation and the MIMO 
relay DCSK-CD system are reviewed and described, respec- 
tively. In Section |nl] performance analysis on the new system 
is carried out, with a close form BER expression derived. 
Simulation results are shown in Section [IV] and conclusions 
are given in Section IVl 

II. System Model 

In this section, the multi-user DCSK system is introduced 
based on Walsh code, and then a new multi-access MIMO 
relay DCSK-CD system is proposed. In the rest of this paper, 
it is assumed that all the independent channels satisfy the 
same channel condition, that is, all channels have the same 
Gaussian-distributed random noise with the power spectral 
density of No. Also, equal-gain combiner (EGC) is adopted at 
the receivers in the systems because it not only can improve 
the error performance but also can be easily implemented |29l . 

A. Principle of DCSK Modulation 

DCSK modulation employs a chaotic signal as the carrier, 
with differential shift keying modulators, for transmission. 
When the chaotic signal is generated by a discrete-time chaotic 
circuit, the global spreading factor is 2/3 ]4J. Fig.[T]shows the 
block diagram of the binary DCSK communication system. 
As seem from this figure, the binary DCSK modulation unit 
transmits a reference segment of the chaotic signal in the 
first half of the symbol duration and, according to the bit 
information, it repeats or reverses the segment in the last 
half of the duration to "1" or "0" respectively. The I th binary 
transmitted signal is represented by a sequence of samples of 
the chaotic signal whose length is 2/3, in which the j th sample 
is given by 

Cj j = 2(1 — l)/3 + 1, • • • , (2/ — l)/3 
(26, - j = (2Z-l)/3 + !,-•• ,2W 



where the I transmitted symbol is denoted by 

bi = {0,1} 



(1) 



(2) 



Here, "0" and "1" appear with the same probability. 

The output signal of the multipath fading channel, repre- 
sented by rj (t), can be expressed as 
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Fig. 1. Block diagram of the DCSK communication system: transmitter (a) 
and receiver (b). 



where on and tj are the independent Nakagami-distributed ran- 
dom variable and the time delay of the i th path, respectively, 
L is the number of fading paths, and rij is the Gaussian- 
distributed random noise with the power spectral density of 
No. The demodulated signal is given by 



= / r(t)r(t-T/2)dt (4) 
It/2 

where T is the duration period of a bit. 



B. DCSK System Based on Walsh Code 

For a multi-access DCSK communication system, the or- 
thogonal Walsh code sequences are adopted to ensure the or- 
thogonality of the signals among different users (3). However, 
the interference among different users can be neglected if the 
Walsh code is synchronized. 

Consider a system with U users. Let / denote the length 
of each carrier segment and 2/3 denote the global spreading 
factor. The 2U -order Walsh code is used to accommodate the 
U users. During the modulation, one keeps / = 2/3/2(7 in 
order to ensure the global spreading factor 2/3 be constant. 
The I th transmitted signal of the u th user can be expressed by 
the 2f/-order orthogonal Walsh code, as follows: 



2U-1 

w 2u -b l ,ic(t 

i=0 



■ T ^ 



< t < T 



(5) 



where bi is the I transmitted symbol and Wk, (k = 
1,2, ••• , 2U) is a row vector of the 2[/-order Walsh code. 
During the demodulation process, the generalized maximum 
likelihood (GML) detector l30l should be applied instead of 
the differentially coherent detector. The block diagram of the 
GML detector for the u th user is plotted in Fig. [2] and the 
weighted energy combined of the / th received signal of the 



user is expressed as 



E, 



u,bi 



T-T/2U 



2U-1 T 

r (*-*--) U; 2«-6 1 ,2(7- l 

. i=0 



dt (6) 



j = 2(1 -l)f3 + !,■■■ ,210 (3) 



Here, r(t — iT/2U) represents the received signal which has 
delay iT/2U. 



Fig. 2. Block diagram of the GML detector for the u th user. 



Fig. 3. The basic model of the 2-user MIMO relay DCSK-CD system. 



For example, to accommodate 2 users, the system should 
use the 4-order Walsh code sequences where 
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(7) 



Each user selects a couple of vectors to form the transmitted 
signal. User 1 can take W\ and u>2 for bit "1" and "0", and 
user 2 can adopt the rest two vectors. 

C Two-user MIMO Relay DCSK-CD System 

Consider a two-hop network system model with two users 
(sources), one destination and one relay. Both the relay and 
destination can employ multiple antennas while each user 
possesses only one antenna. The channels constructed by the 
transmit-receive antenna pairs are independent and subjected 
to static block frequency-selective fading, meaning that the 
channel state remains constant during each transmission pe- 
riod. In this system, the users do not cooperate with each other 
in sending messages but only the relay helps them to transmit 
messages to the destination. To simplify the analysis, it is 
assumed that a transmission period is divided into broadcast 
phase and cooperative phase, denoted as 1 st time slot and 2 nd 
time slot, respectively. In the 1 st time slot, users broadcast 
their messages to other terminals (relay and destination), and 
then the relay cooperates with two users to send their messages 
to the destination in the 2 nd time slot. The basic system model 
is presented in Fig. [5] where each user can support only 
one transmit antenna, the relay possesses M R antennas used 
under ideal condition (error free (EF) at the relay antennas) 
or decode-and-forward (DF) protocol [371], and the destination 
has Mo receive antennas. The transmission principle of the 
system is shown in Fig.@] where X RA (t) and X RB (t) repre- 
sent the information reconstructed by the relay. In particular, 
X RA (t) = X A (t) and X RB {t) = X B (t) for the EF situation; 
for the DF case, X RA (t) = X A (t) and X RB {t) = X B (t) if 
the relay is able to decode the signal correctly, and X RA (t) = 
X RB {t) = otherwise. The function F(X RA (t),X RB (t)) 
describes the processing of the reconstructed signals. In the 
proposed system, F(X RA (t),X RB (t)) = X RA (t) + X RB (t). 
Consequently, the received signals at the relay and destination 
antennas can be expressed as follows: 
1 st time slot: 

r R (t) = H AR ® X A (t) + H BR ® X B (t) + Z R (t) 
r D (t) = H AD <8> X A {t) + H BD <g> X B (t) + Z D {t) 
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Fig. 4. The transmission principle of the proposed system. 



2 nd time slot: 

Hrd ( 



"?>(*) 



(X RA (t)+X RB (t)) + Z 2 D (t) 



(9) 



Here, the superscripts "1" and "2" are used to denote the 1 st 
time slot and the 2 nd time slot, <g) represents the convolution 
operator, r D (t) and r 2 D (t) are the total received signals at 
destination for users A and B in the 1 st time slot and the 
2 nd time slot, respectively, r R (t) is the total received signals at 
relay for two users in the 1 st time slot. X A (t) and X B (t) is the 
signal transmitted by user A and B, Z R (t) and Z D (t) (i = 1, 2) 
are white Gaussian noise random variables with zero mean and 
the power spectral density Nq. The channel multipath impulse 
response H mn are modeled by a linear time-invariant process, 
i.e., H rnn = Y^=i a i${t — T i), where on and r, denote the 
attenuation and time delay of the i th path, respectively, and L 
represents the number of the multipath taps. 

The channel coefficients | on \ follows a Nakagami distri- 
bution, with the probability density function (PDF) given as 



/|ai|( C 



r(mj) 



a 



2m; -1 



cxp 



(10) 



where m; is the fading parameter of the i -path, Qi equals 
E(af), and T(-) denotes the Gamma function. To simplify the 
analysis, it is assumed that each fading channel has a uniform 
scale parameter, meaning that Vli/mi is constant for all the 
paths. Meanwhile, the receiver knows all multipath parameters, 

thus Ef=i^i = Ef=i^(«?) = i. 

As two conventional DCSK systems, the DCSK-NC system 
does not use any cooperative technique; for the DCSK-CC 
system fPHl . which involves no relays, the users do cooperate 
with each other to transmit messages. Fig. [5] shows the basic 
model of the 2-user DCSK-CC system, where the parameters 
are shown in the first paragraph of this subsection. Referring 
to this figure, one can observe that the DCSK-CD system (see 
Fig. [3} is quite different from the DCSK-CC system. In the 
DCSK-CC system, the users broadcast their messages to other 
terminals (the other user and destination) in the 1 st time slot. 
Then, these two users cooperate with each other to send their 
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Fig. 5. The basic model of the 2-user DCSK-CC system. 



messages to the destination through two different channels in 
the 2 nd time slot. 

Comparing with the DCSK-CC system, a pair of transmit- 
receive antennas can be used to construct an independent 
fading channel in the proposed DCSK-CD system without 
increasing much complexity. Thus, the anti-fading capacity 
can be significantly enhanced. As a result, a considerable 
performance gain is obtained by incorporating the MIMO relay 
technology into the DCSK system. 

III. Performance Analysis 

In this section, the exact BER expression of DCSK-NC 
system l32l is re-derived over Nakagami fading channels in 
order to obtain the bit error rate (BER) formula about the 
multi-access MIMO relay DCSK-CD system. Furthermore, a 
method is presented in closed-form for approximating the BER 
expressions of the DCSK-NC and the proposed system. In 
the following analysis, we assume that the energy per bit is 
constant, denoted by E b , the Walsh code used in the system 
has been perfectly synchronized, and the time delay of each 
path Ti is much shorter than the bit duration (r^ <C 2f3T s , T s 
is the sampling period of the chaotic signal), so that the inter- 
symbol interference (ISI) in the transmission can be neglected. 

A. Exact BER of Conventional DCSK-NC System 

First, the BER expression of the DCSK-NC system is re- 
derived over Nakagami fading channels. 

The signal-to-noise ratio (SNR) per bit of the j b at the 
receiver can be expressed as 

L L 



lb 



(11) 



where the bit energy Eb is defined by Et, = 2j3E[(cj) 2 ], 
Nq is the power specral density of noise. L is the number 
of independent fading paths, ji is the instantaneous SNR of 
the i th path, and m is the Nakagami-distributed varaible with 
parameter rm. 

It is well known that the square of a Nakagami-distributed 
random variable follows a gamma distribution. The PDF of a 
gamma-distributed variable X with parameters a > 0, b > 0, 
denoted as X ~ G(a, b), is given by 

a—l —x/b 

/(*) = ^p *>0 d2) 

As ctj is the Nakagami-distributed variable with parameter m^, 
we have 7i = (E b /N )a% ~ G (rm, (E b /N )Sli/mi)- The 



MGF of the variable X ~ G(a, b) is given by ED 

M(s) = (1 - sb)- a (13) 

Since 71, • • • , 7 l are statistically independent and f2j/mj is 
constant, we can get the MGF of j b = Ylf=i 7» through 

(14) 



M iA s ) = H(i - s — ) 



Thus the distribution of -f b can be obtained using the inverse 
MGF transform of ([Pil l, resulting in 



lb-G 



E b J2^ 



\ 



m. 



G mL 



E b /N 
mL 



(15) 
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where mi(m\ = m,2,- ■■ ,= = rn) is the Nakagami 
fading parameter of the i th path which controls the depth of 
the fading envelope. Hence, the PDF of j b can be written as 



f{lb) 



lb 



mL-1 Tb/^ mL ) 



( E b /Ng 

\ mL 



mL 



(16) 



T{mL) 



The conditional BER of the DCSK-NC system as a function 
of the received SNR is given by l32l 



BER( 76 ) = i erfc 



4 

76 



1 



J_ 

27fc 




(17) 



76 



27fc + / 



where / is the length of each carrier segment, erfc(-) is defined 
as 



erfc(x) 



"dA 



(18) 



and Q(-) is the Gaussian Q-fuction. Combining (TToT l and 
dTTI i, the BER of the DCSK-NC system can be obtained by 
averaging the conditional BER, as follows: 



BER = / BER( 7h )/( 76 )d7 6 
'0 



(19) 



This formula will be used for deriving the exact BER expres- 
sions of the MIMO relay DCSK-CD system with different 
relay protocols in the next section. 

B. Exact BER of Multi-access MIMO Relay DCSK-CD System 

In the multi-access MIMO relay DCSK-CD system, assume 
that the fading coefficient on of each path is independent and 
identically-distributed (i. i. d.). Consider a system includes n 
users with only one antenna, one receiver with Md antennas 
and one relay with Mr antennas. To simplify the analysis, also 
assume that the transmitted messages arrived at each receive 
antenna possess identical energy and all independent channels 
are subjected to the same channel condition. Moreover, let 2/3 
denote the global spreading factor, and the dsD, dsn, d^D 
denote the distance of source to destination (S-D), source to 
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relay (S-R) and relay to destination (R-D) links, respectively, 
thus <1sd '■ dsR '■ <Ird represents the geometric positions 
of the users. The path loss of the S-D link is given by 
PL(S,D) = l/d 2 SD . 

For easily deriving the the exact BER of the multi-access 
MIMO relay DCSK-CD system with different protocols, we 
should introduce the following two theorems. 

Theorem 1. Given N independent gamma random variables 
X U X 2 ,-- ,X N , where X k ~G(a k ,b) (k = l,-- ,7V). The 
sum of these variables X = 2~2^=i also follows a gamma 
distribution, denoted as X ~ G(^2 k=1 a k , b). 

Theorem 2. Assume that N independent gamma ran- 
dom variables X\,X 2 , — - ,Xn are distributed as X k ~ 
G(a k ,b k ){k = ,N) where b x ± b 2 ,-- b N . The 

PDF of the sum of these variables X = 2~2 k =i ^ k * s ex P resse d 
as 



shown as 



oo 
i=0 

which is subjected to 

N 



(20) 



c = i[(bo/b k r* 

k=l 

Vt+i = — — r y~)tztT}i+i-t i = 0,1,2, •• 
t=l 

N 

Zj =J2ak(l-b /b k y/j j = l,2,-- 



(21) 



k=l 
N 

p = ^2a k >0 

k=l 

b () = mm(b k ) 

k 



where 770 = 1. 

1 ) BER with EF protocol for the relay ( ideal condition ): 
Assume that the relay can receive messages from the users 
without errors, namely, the BER of S-R link is zero. Obviously, 
there are Mo, Mr and Mr Mo independent Nakagami fading 
channels for S-D link, S-R link and R-D link. The conditional 
BER and total BER expressions of the proposed system are 
expressed by $1% and $1% . Thus, what needs to do is to 
deduce the PDF of the total received SNR. Let the received 
SNRs for the relay and the destination in the 1 st time slot be 
denoted by 7,5 ^ and jsd, and similarly the received SNR for 
the destination in the 2 nd time slot be denoted by Jrd- 

In the 1 st time slot, namely the broadcast phase, based 
on the equivalent transmitted SNR of this channel is 
(Eb/2Mo)/No, one can easily obtain the distribution of the 
received SNR for the fc th source-destination antenna pair, as 
follows: 

ISD-k - G (ml, (E b /N Q )/(2M D d 2 S omnL)) 
Thus, the distribution of jsd can be deduced using TheoremQ] 



Mr 



Isd = ^2 ISD-k ~ G I M D mL, 



k=l 



E b /N 



2M D d 2 SD mnL J (22) 



G(aiA) 



where m is the Nakagami fading parameter of each 
path, n is the number of users. In the following, a\ 
and 61 are used as short-hand notations for MomL and 
(E b /N )/(2M od 2 SD mnL), respectively. 

Similarly, in the 2 nd time slot, one can get the distribution 
of received SNR for the k th antenna pair for R-D link, 
namely ^ RD -k ~ G (mL, (E b /N )/(2M R Mod RD mnL)). 
Thus, ^Ro is given by 



M R M D 

1RD = ^ jRD-k 
k=l 



G I M R M D mL 



G(a 2 ,b 2 ) 



E b /N 



2MRMod 2 RD mnL 



(23) 

Here, likewise, a 2 and b 2 are short-hand notations for 
M R M D mL and {E b /N ) / 1 (2M R M D d 2 R D mnL). 

Finally, using the equal-gain combining method, the total 
received SNR for destination with two signals from the relay 
and one source jo can be formulated as 



Id = (isd + 1rd)/V2 



(24) 



where dividing by s/2 (normalized factor) is used to keep the 
total received energy for the destination from two different 
terminals in the proposed system be consistent with the ones 
of the DCSK-NC system. 

Combining (f22b-(f24b with Theorems [T] and [2] we can write 
the PDF of jo as 



firm) = 



^,a 1 +a 2 -l e - 7D /(b 1 /V2) 

(6 1 /V2) ai + a3 r( ai +a 2 ) 



i=0 



T(p + i)b p 



P+i 



if 61 = b 2 



if 61 ^ b 2 



(25) 



where all parameters are as shown in (T2TI)-(l23l). 

Substituting (T% and (l25l l into ( fT~9b , one can get the exact 
BER formula of the proposed system under the ideal condition 
as d26l i. where / is the length of each carrier segment. 

2) BER with DF protocol for the relay: As one of the 
most popular relaying techniques employed in cooperative 
diversity systems, DF protocol performs such that if the relay 
can decode the received signal correctly, then it forwards 
the decoded message to the destination in the 2 nd time slot, 
otherwise, it remains idle. The BER of the DF case is based 
on the one of the EF protocol, described as 

BERdf = BERg^ • BERg^i + (1 — BER5/;) ■ BER^) (27) 

where BER Sfl , BER 5 d and BER^ denote the BER of the 
relay receiver with one signal from one source, the BER of the 
destination receiver with one signal from one source, and the 
total BER of the destination receiver with two signals from 
one source and the relay, respectively. The distribution of the 
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BERwp — BERn — 



= < 



BER( 7D )/( 7D )d 7D 

oo ^i+a 2 -l e - 7D /(fc 1 / v / 2) 

o (6 1 /V2) 0l+02 r(a 1 + a 2 ) "* \ V 2 lD + f 



oo / 00/ n-f-i— 1 — Tr)/bn 



i=0 



p+i 




d 7 £) if 61 = b 2 



djo if 61 7^ b 2 



{b 1 /V2) ai+a2 T(a 1 +a 2 ) J 



C 



E 



7 ai+a 2 -lg-7D/(6i/V2) 




d 7 c if &i = 62 



d 7 D if &i ^ &2 



(26) 



received SNR for the relay is 



1SR = ^2 "fsn-k ~ G I M R mL, 



fc=l 



E b /N 



2M R d 2 SR mnL 



(28) 



= G(a 3 ,6 3 ) 



where jsR-k is the received SNR for the fc th S-R channel, 
expressed as -f S R-k ~ G (mL,(E b /N )/(2M R d 2 SR mnL)), 
and and 03 and 63 are short-hand notations for M R mL and 
{E b /N )/(2M R d 2 SR mnL). 

The distributions of 7 5d and 7 £> here are the same as the 
ones under the ideal condition, with expressions shown in (l22l 
and (l25T l respectively. Then we can obtain the PDFs of these 
three variables, denoted by /( 7 s_r)> /( 7 sd), and /( 7 d). The 
conditional BERs for these variables are identical, as shown in 
(fTTT i. Accordingly, the expression of BER/3 is shown in d26i i, 
while BERgfl and BERg^ can be calculated by ( fT9] l, thus 
we have 



BERgij 



00 7|£V 



&3 3 r(a 3 ) y 



a 3 — l„ — JSR/b3 
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Finally, substituting <|26), (|29j and ([30} into (|27j, the BER of 
the proposed system with DF relay protocol can be formulated. 

C. Approximate Results of the DCSK-NC System and the 
Proposed System 

Although the exact average BER can be computed through 
numerical evaluation of the double integral in (fT9l with 
different PDFs of received SNRs. It is very difficult to 



calculate the integrations of ( fT9l and ( T26l l arising from the 
nonlinear behavior of the function. To resolve this problem, 
an approximated method based on the MGF l33l is proposed 
to simplify the expression of the BER into a closed-form. 

1 ) DCSK-NC system: To begin with, define w as a function 
of jb, as 

w = (31) 
2 7 fc + / 

Substituting OTb into $1% gives the BER expression via 
averaging the conditional BER of w 
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where the Gaussian Q-function may be represented as 
Q(V^>) = (l/n)jf 2 e-^ed6. 

Since the MGF of w is the Laplace transform of f w (w) 
with the exponent reversed in sign, defined as M w (s) = 
j°o e sw f w ( w }d W} putting it into d32l yields 

r ir/2 
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(33) 



Assume that -f b ~ G(a, b). Then, the variable w can be ap- 
proximately described by a Gamma distribution 1351 , denoted 
as G(a w , b w ), with parameters given by 

a w = a _ 

V7+/; (34) 

&7» = 



2( 7 + //2)3 

where 7 = £f =1 7i - (^b/^o) E.ti = E b /N and 

/ is the length of each carrier segment. 

The MGF of w ~ G(a w , b w ) can be formulated using ( TT3b 

M w («) = (1 - S 6„,)" atu (35) 
Substituting d35l l into d33l gives the following expression: 
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which can be written in closed-form by using the definite 
integral derived in 1331 . Hence, the approximate result of the 
BER formula of the DCSK-NC system is 
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(37) 



where 2-fi("j S S •) is the Gaussian hypergeometric function 
shown in (36). It should be noted that the approximate result 
will become more accurate as the number of paths increases. 

2) MIMO DCSK-CD system: We only consider the approxi- 
mate BER expression with EF protocol over the Nakagami fad- 
ing channel, which is the foundation to derive the one with the 
DF case. It is known from d25l l that jjj = (jsd + 7hd)/v / 2 
is a gamma-distributed variable if bi = b^, denoted as 
7d ~ G(ai + a,2,bi/\/2). Combining the equation bi = 62 
with ( 122b and (1231) . we conclude that 70 follows gamma 
distribution provided that dsD/duD = V^V/r- In this case, 
one can also use (f37b to estimate the BER of such a system 
with different parameters (a e f,6 e f) where 

2 



(ai + a 2 ) 



flD + f/2 



b c f = 



V 1D + f 

bi^pjjp + /) 2 
2%/2( 7 d + //2) 3 



(38) 



Here, ~/ D = J2 l=1 lw = E(*y D ). 

Hence, the approximate BER expression is 
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(39) 



2 + b r 



The approximate result of the DCSK-CD system with DF 
protocol can be obtained similarly. However, in the case of 
bi 7^ bi, it is not possible to derive the MGF for since 
it does not follow gamma or any other regular distributions. 
Therefore, the close-form BER expression is unavailable to- 
day. 

IV. Simulation Results 

In this section, some numerical calculation, approximation 
and simulation results of BER are presented. Since the perfor- 
mance of the multi-access system with more than three users 
is the same as that with two users when the energy of each 
user is constant, we only consider the systems with two users 
here. 

It is assumed that all channels have the same Gaussian- 
distributed random noise and the same path-loss model. The 
transmission energy per bit is kept constant and the energy of 
each path is allocated uniformly. Unless otherwise stated, the 
simulations are performed in 2-user systems over Nakagami 
fading channels, with parameters m = 1 and L = 2. The 
delay vector t = (ti,t-2) = (0,T S ). The global spreading 
factor of DCSK is much higher than the path delay r^, 
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Fig. 6. BER curves of three different systems with the EF protocol (a) and 
the DF protocol (b). 



namely 2/3T s ^ T{. In the experiments, (Mr,Md) is used 
to represent the number of antenna pairs at the relay and the 
destination, respectively. 

A. Performance Comparison among the MIMO Relay DCSK- 
CD, DCSK-CC and DCSK-NC Systems 

Fig. |6] plots the BER curves of the DCSK-NC, the con- 
ventional DCSK-CC and the MIMO relay DCSK-CD systems 
with EF (the ideal condition) protocol and DF protocol over 
Nakagami fading channels, where d$D ■ dsn ■ d^o and 
(Mr,Md) are set to 1 : 1 : 1 and (1,1), respectively. The 
global spreading factor considered here is 2/3 = 128 (/ = 32). 
Referring to this figure, the BER performance of the proposed 
DCSK-CD and the DCSK-CC systems is remarkably better 
than the DCSK-NC system at a BER of 4 x 1CT 3 . It also shows 
that the performance of the DCSK-CD system outperforms 
the DCSK-CC system by about 1 dB for BER = 4 x 1CT 3 
under these two conditions. Furthermore, we observe that 
the numerical BER and the simulation results are highly 
consistent. More performance gain of the proposed system can 
be expected as (Mr,Md) becomes larger. 
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Fig. 8. BER curves of the MIMO relay DCSK-CD systems with different pig. 10. BER curves of the proposed systems with different spreading factors, 
antenna pairs. 



Further, Fig. [7] plots the BER curves of the proposed 
DCSK-CD system with EF and DF protocols for different 
fading depths. The parameters used are the same as those 
in Fig. [6] except that the fading parameter m now varies 
(m = 0.5,0.8,1,2). Referring to this figure, the performance 
gap (distance between the two curves for a given BER) 
between the EF and DF protocols becomes larger as the fading 
depth increases (m decreases), for example, the EF protocol 
has a performance gain about 1.5 dB over the DF case at a 
BER of 2 x 10~ 3 when m = 0.5. However, we also observe 
that the performance of the proposed system with DF protocol 
very approaches that of the EF case when the value of m is 
over one because the relay can successfully decode most of 
the transmitted bits. 

B. Performance of Two-user MIMO Relay DCSK-CD System 

To further justify the applicability of the proposed system, 
some practical environments are now discussed for this system 
with DF protocol. Assume that dsD ■ dsR '■ dun = 1:1:1 
and 2/3 = 128. The BER curves for different (Mr, Md) 
pairs are shown in Fig. [8] One can observe that the BER 



performance of the proposed system will be improved as 
the diversity order MrMd increases for obtaining more 
spatial diversity gain. For example, the proposed system with 
(Mr,Md) = (3,2) has a gain about 3 dB over that with 
(1,1) at a BER of 2 x 10~ 4 , and it is likely to be even more 
in the higher SNR region. Also, the influence of different 
distance ratios on the system performance was investigated. 
The users always choose a relay which is much closer to the 
destination and cooperate with them to transmit messages. In 
this situation, we assume that 2/3 and (Mr, Md) here maintain 
128 and (2,2). The distance ratio dsD ■ dsn. '■ <1rd is set to 
be 1 : 0.8 : 0.4 and 1:1:1, respectively. The BER curves of 
these two cases with DF are plotted in Fig. [9] It illustrates that 
when dsD is fixed, the system with less dsn or dan distance 
exhibits better performance for obtaining extra path-gain (less 
path-loss). 

In addition, Fig. [10] depicts the BER curves of the systems 
with three different spreading factors. The parameters dsD ■ 
dsn '■ dRD, (Mr, Md) and m are assumed to be 1 : 1 : 
1, (2, 2) and 1 while the global spreading factors 2/3 = 64, 128 
and 256 are considered. In this figure, we observe that better 
performance can be expected for a smaller spreading factor 
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when the energy per bit is constant provided that 2/3T s ^> n 
since the average energy of each sampling point E{(ci) 2 } = 
Eb/{20) will decrease as 2/3 becomes larger and thus results 
in worse anti-noise capacity. 

C. Discussion about the Threshold of Antenna Numbers for 
the Proposed System 

Given a certain value of fading factor m and path number 
L of the Nakagami fading channel, the diversity order of the 
proposed system is MrMd- As the energy per bit is set to 
constant, i.e., Ef, is unit energy, the transmission energy of 
each individual channel constructed by transmit-receive an- 
tenna in the R-D link is Ei,/(2MrMd) (another Eb/2 energy 
is used by source in the 1 st time slot). Consequently, the 
system performance should be improved as MrMd increases 
when this parameter is small since the performance gain 
from a higher diversity order is larger than the performance 
loss from the worse anti-noise capacity, especially in the 
high SNR region. However, the system performance is no 
longer improved if MrMd exceeds a certain threshold as the 
transmission energy is extremely low so that the effect of the 
diversity gain can just offset the effect of the deteriorated anti- 
noise capacity. Consequently, there is a threshold of the value 
of MftMrj in the proposed system where the influence of 
parameter Mr is equivalent to that of Mr>, and this threshold 
can be obtained through observation of the corresponding 
theoretical performance curves. 

Based on the aforementioned analysis, the system perfor- 
mance will become better as the diversity order MrMd be- 
comes larger, but the complexity will increase simultaneously. 
Moreover, the performance gain is negligible when MrMd 
exceeds a certain threshold. So, a trade-off between system 
complexity and performance should be considered. As the 
effect of the parameters Mr and Mrj are equivalent in such a 
system, the discussion here is only focusing on the influence 
of the parameter Mr while keeping Mrj constant, and vice 
versa. Assume that dsD ■ d$R. : <Ird = 1:1:1, 2/3 = 128 
and Mrj = 2. Fig. Qj] shows the exact BER curves obtained 
from numerical calculation with DF protocol in Sect. IIII-BI 
Referring to this figure, the system performance is almost 
no longer improved if the Mr is over four, because the 
transmission energy for the channel that constructed by the 
transmit-receive antenna pair is extremely low, which makes 
the performance gain from the spatial diversity merely the 
same as the performance loss from the worse anti-noise 
capacity. Accordingly, M R = 4 and (MrMd)* = 8, can 
be seen as a threshold for Mr and MrMd in the proposed 
system, respectively. 

D. Performance Comparison between Approximate and Exact 
BER Results 

From the numerical and simulation results shown in the 
previous subsections, we conclude that they agree quite well. 
In this subsection, the performance of the approximate (APP) 
BER expressions in Sect. IIII-CI is compared with the exact 
BER results in Sect. IIII-AI and IIII-BI thereby verifying the 
accuracy of the closed-form expressions. 
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Fig. 1 1 . The BER curves of the proposed system with different numbers of 
relay antennas. The number of destination antennas (Afo) is set to 2. 

We assume that the system parameters d$D '■ d$R ■ 
<1rd, (M R ,M D ) and 2/3 are set to 1 : 1 : 1, (1, 1) and 128, 
respectively. Therefore, the precondition of the formula ( |39l ) 
is fulfilled, i.e., d SD /d RD = v/Mr = 1. The APP BER and 
the exact BER curves of DCSK-NC and MIMO relay DCSK- 
CD systems with different numbers of paths are presented 
in Fig. [12] It implies that as the number of paths increases, 
the performance curves of approximate BER and exact BER 
expressions become more consistent for such two systems. In 
particular, referring to Fig.Q~2jb), the APP BER curve is much 
closer to the exact one when L = 8 as compared to those of 
L = 2 and L = 4. Based on the experimental results, we may 
conclude that the approximate BER expressions are expected 
to be more accurate with a large number of paths. 

V. Conclusions 

In this paper, a MIMO relay DCSK-CD system has 
been proposed as a comprehensive cooperation scheme to 
strengthen the robustness of anti-fading in a wireless network, 
where both the relay and the destination can support multiple 
antennas, to get more spatial diversity gain. An exact BER 
formula and the corresponding closed-form approximate ex- 
pression are derived over Nakagami fading channels, which 
agree well with the simulation results. Comparing to the con- 
ventional DCSK-NC and DCSK-CC systems, the performance 
gain of the proposed system with one antenna for relay and 
destination is about 4 dB and 1 dB at a BER of 4 x 1CT 3 , 
respectively, and more performance gain can be expected with 
a higher diversity order MrMd- Moreover, it was observed 
that a trade-off between system performance and complexity 
is needed when the energy of each user is kept constant, and 
an optimal value of the number of antennas at relay (M R = 4) 
is obtained through observation of the theoretical performance 
curves when Mo is set to 2. According to these advantages, 
the proposed system has great potential in improving the 
performance of energy-constrained wireless networks, such 
as low-power and lower-cost WPANs. In the future, we will 
strive to derive a close-form BER formula for the case of 
dsD /dRD 7^ VMr thereby resolving the remainder practical 
problem in introducing the proposed system. 
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Fig. 12. The APP and the exact BER curves of the DCSK-NC system (a) 
and the MIMO relay DCSK-CD system (b) with different numbers of paths. 
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